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Abstract 

Given a set S of spheres in E d , with d > 3 and d odd, having a fixed number of m distinct 
radii p±,p2,.-. ,pmi we show that the worst-case combinatorial complexity of the convex hull 

— 

CHd{^) of £ is ®(^2x<ijtj<m n i n j 2 )) where n, is the number of spheres in £ with radius 
Pi. Our bound refines the worst-case upper and lower bounds on the worst-case combinatorial 
complexity of Ci?d(£) for all odd d > 3. 

To prove the lower bound, we construct a set of 0(ni + n.2) spheres in E d , with d > 3 
odd, where n» spheres have radius pi, i = 1,2, and pi ^ pi, such that their convex hull has 

I — I I — I 

combinatorial complexity Q(nin 2 2 + n2 r n 1 2 ). Our construction is then generalized to the case 
where the spheres have m > 3 distinct radii. 

For the upper bound, we reduce the sphere convex hull problem to the problem of computing 
the worst-case combinatorial complexity of the convex hull of a set of m d-dimensional convex 
polytopes lying on m parallel hyperplanes in ~E d+1 , where d > 3 odd, a problem which is of 
independent interest. More precisely, we show that the worst-case combinatorial complexity 

of the convex hull of a set {Pi, 7^2, • • • ,P m } of m d-dimensional convex polytopes lying on m 

1 — 

parallel hyperplanes of E +1 is 0(J2i<i^£j< m n i n j )> where rii is the number of vertices of Vi- 
This bound is an improvement over the worst-case bound on the combinatorial complexity of 
the convex hull of a point set where we impose no restriction on the points' configuration; using 
the lower bound construction for the sphere convex hull problem, it is also shown to be tight 
for all odd d > 3. 

Finally: (1) we briefly discuss how to compute convex hulls of spheres with a fixed number 
of distinct radii, or convex hulls of a fixed number of polytopes lying on parallel hyperplanes; 
(2) we show how our tight bounds for the parallel polytope convex hull problem, yield tight 
bounds on the combinatorial complexity of the Minkowski sum of two convex polytopes in E d ; 
and (3) we state some open problems and directions for future work. 

Key words: high-dimensional geometry, discrete geometry, combinatorial geometry, combinato- 
rial complexity, convex hull, Minkowski sum, spheres, convex polytopes, parallel hyperplanes 
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1 Introduction and results 



Let £ be a set of n spheres in K d , d > 2, where the dimension d is fixed. We call II a supporting 
hyperplane of £ if it has non-empty intersection with £ and £ is contained in one of the two closed 
halfspaces bounded by II. We call H a supporting half space of the set £ if it contains all spheres in £ 
and is limited by a supporting hyperplane IT of £. The intersection of all supporting halfspaces of £ 
is called the convex hull Ci?d(£) of £. The definition of convex hulls detailed above is applicable not 
only to spheres, but also to any finite set of compact geometric objects in K d . In the case of points, 
i.e., if we have a set P of n points in K d , the worst-case combinatorial complexitj0 of CHd(P) is 
known to be Q(n^^). Moreover, there exist worst-case optimal algorithms for constructing CHd(P) 
that run in 0(n^i + nlogn) time, e.g., see [T6l l25l l2l IM 113]. Since the complexity of CHd(P) 
may vary from 0(1) to 0(nW), a lot of work has been devoted to the design of output-sensitive 
algorithms for constructing CHd(P), i.e., algorithms the running time of which depends on the size 
of the output convex hull CH d (P), e.g., see {L2\ (H [271 US Ell E2 CD EJ \W\ EE] ■ For a nice overview 
of the various algorithms for computing the convex hull of points sets, the interested reader may 
refer to the paper by Erickson |15| . while Avis, Bremner and Seidel [3] have a very nice discussion 
about the effectiveness of output-sensitive convex hull algorithms for point sets. 

Results about the convex hull of non-linear objects are very limited. Aurenhammer [3J showed 
that the worst-case complexity of the power diagram of a set of n spheres in E^, d > 2, is 0(n^^), 
which also implies the same upper bound for the worst-case complexity for the convex hull of the 
sphere set. Rappaport |26| devised an O (nlogn) algorithm for computing the convex hull of a set 
of discs on the plane, which is worst-case optimal. Boissonnat et al. [6] give an 0(n^l + nlogn) 
algorithm for computing the convex hull of a set of n spheres in K d , d > 2, which is worst-case 
optimal in three and even dimensions, since they also show that the worst-case complexity of the 
convex hull of n spheres in E 3 is 0(n 2 ). Finally, their results hold true for the case of homothetic 
convex objects. Boissonnat and Karavelas [7J settled a conjecture in [6]: they proved that the worst- 
case complexity of a set of n spheres in K d , d > 2, is 0(n^), which also implied that the algorithm 
presented in [6] is optimal for all d. As far as output-sensitive algorithms are concerned, Boissonnat, 
Cerezo and Duquesne [5] showed how to construct the convex hull of a set of n three-dimensional 
spheres in 0(nf) time, where / is the size of the output convex hull, while Nielsen and Yvinec |23] 
discuss optimal or almost optimal output-sensitive convex hull algorithms for planar convex objects. 

In this paper we consider the problem of computing the complexity of the convex hull of a set 
of spheres, when the spheres have a fixed number of distinct radii. This problem has been posed by 
Boissonnat and Karavelas [7J, and it is meaningful for odd dimensions only: in even dimensions the 
complexity of both the convex hull of n points and the convex hull of n spheres is 0(n^2 J ) = G(n^T ) ; 
i.e., the two bounds match. 

Consider a set of n spheres £ in K d , where d > 3 and d odd, such that the spheres in £ have a 
fixed number m of distinct radii p\,p2, ■ ■ ■ , p m . Let n, be the number of spheres in £ with radius 
Pi. We say that p\ dominates £ if n\ = Q(n). We further say that £ is uniquely (resp., strongly) 
dominated, if, for some A, p\ dominates £, and nj = o(n) (resp., Hi = 0(1)), for all i 7^ A. Using this 
terminology, we can qualitatively express our results as follows. Firstly, if £ is strongly dominated, 
then, from the combinatorial complexity point of view, CHd(T,) behaves as if we had a set of points, 
or equivalently a set of spheres with the same radius. If, however, £ is dominated by at least 
two radii, CH^iT,) behaves as in the generic case, where we impose no restriction on the number of 
distinct radii in £. Finally, if £ is uniquely dominated (but not strongly dominated), the complexity 

In the rest of the paper, and unless otherwise stated, we use the term "complexity" to refer to "combinatorial 
complexity". 
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of C-ffrf(S) stands in-between the two extremes above: the complexity of C-ffd(E) is asymptotically 

larger than the points' case (or when we have spheres with the same radius), and asymptotically 

smaller than the generic case, where we impose no restriction on the number of distinct radii in 

E. From the quantitative point of view, our results refine the results in [7] for any odd dimension 

I i\ 

d > 3 as follows. We prove that the worst-case complexity of CHd{S) is 6(Ei<j^j<m n * n j 2 )■ This 
tight bound constitutes an improvement over the generic worst-case complexity of C-ff^E) if E is 
uniquely dominated, since in this case the complexity of CH^T,) is o(n^z^). On the other hard, it 
matches the generic worst-case complexity of C-ffd(E), if E is dominated by at least two radii: in this 
case the worst-case complexity of CHd(Y,) becomes ©(n^l). Finally, if E is strongly dominated, 
the complexity of CH^T,) is 0(71^2-1), i.e., it matches the worst-case complexity of convex hulls of 
point sets (or sets of spheres where all spheres have the same radius). 

To establish the lower bound for the complexity of C-ffd(E), we construct a set E of 0(ni + 77,2) 
spheres in K d , for any odd d > 3, where n\ spheres have radius p\ and ri2 spheres have radius 

L-J L-J 

P2 7^ Pi, such that worst-case complexity of CHd(T,) is Q(nin 2 2 + ). This construction is 

then generalized to sets of spheres having a fixed number of m > 3 distinct radii. More precisely, 
we construct a set E of n = Y^=i n i spheres, where ra^ spheres have radius pi, with the p^s being 

L-J 

pairwise distinct, such that the worst-case complexity of Ci^rf(E) is ^(Xa<i^/<m n i n j 2 )• 

To prove our upper bound we use a lifting map, introduced in |6], that lifts spheres Uj = (q, rj) in 
K d to points pi = (q, rj) in K d+1 . The convex hull Ci^(E) is then the intersection of the hyperplane 
{xd+i = 0} with the Minkowski sum of the convex hull CHd+i{P) and the hypercone Ao, where 
P is the point set {pi,P2, ■ ■ ■ ,Pn} in IE d+1 , and Ao is the lower half hypercone with arbitrary apex, 
vertical axis and angle at the apex equal to j. When the spheres in E have a fixed number m of 
distinct radii, the points of P lie on m hyperplanes parallel to the hyperplane {x,i+i =0}. In this 
setting, computing the complexity of Cifd(E) amounts to computing the complexity of the convex 
hull of m convex polytopes lying on m parallel hyperplanes of K d+1 . This observation gives rise 
to the second major result in this paper, which is of independent interest, and gives as corollary a 
tight bound on the worst-case complexity of the Minkowski sum of two convex d-polytopes. Given 
a set = {V\,V2, ■ ■ ■ ,Pm} of m convex <i-polytopes in E d+1 , with d > 3 and d odd, we show 

I -J 

that the worst-case complexity of the convex hull CHd+i(& > ) is 0(X^i<i^<m n i n j 2 )> where raj is 
the number of vertices of V%- Our upper bound proof is by induction on the number m of parallel 
hyperplanes. The lower bound follows from the lower bound on the complexity of the convex hull 
of spheres having m distinct radii. Our bound constitutes an improvement over the worst-case 
complexity of convex hulls of points sets, if a single polytope of 8? has 0(ra) vertices, whereas all 
other polytopes have o(ra) vertices, where n is the total number of vertices of all m polytopes, while 
it matches the worst-case complexity of convex hulls of points sets if at least two polytopes have 
0(n) vertices. 

The rest of our paper is structured as follows: In Section [2] we detail our inductive proof of the 
upper bound on the worst-case complexity of the convex hull of convex polytopes lying on parallel 
hyperplanes, while in Section [3] we discuss how to compute this convex hull. In Section U we prove 
our upper bound on the worst-case complexity of the convex hull of a set of spheres. Next we present 
our lower bound construction for any odd d > 3 in two steps: first for sphere sets with two distinct 
radii and then for sphere sets with m > 3 distinct radii. We end the section by discussing how this 
lower bound yields a tight lower bound for the problem of the previous section. In Section [5] we 
explain how to modify the algorithm by Boissonnat et al. [6] so as to almost optimally compute 
the convex hull of a set of spheres with a fixed number of distinct radii. Finally, in Section [6] we 
summarize our results, we explain how our results yield tight bounds for the complexity of the 
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Minkowski sum of two convex polytopes, and state some open problems. 

2 Convex hulls of convex polytopes lying on parallel hyperplanes 

A convex polytope, or simply polytope, V in K d is the convex hull of a finite set of points P in K d . A 
polytope V can equivalently be described as the intersection of all the closed halfspaces containing 
P. A face of V is an intersection of V with hyperplanes for which the polytope is contained in one 
of the two closed halfspaces determined by the hyperplane. The dimension of a face of V is the 
dimension of its affine hull. A fc-face of V is a /c-dimensional face of V . We consider the polytope 
itself as a trivial ci-dimensional face; all the other faces are called proper faces. We will use the term 
d-polytope to refer to a polytope the trivial face of which is ti-dimensional. For a (i-polytope V, the 
0-faces of V are its vertices, the 1-faces of V are its edges, the (d — 2)-faces of V are called ridges, 
while the (d— l)-faces are called facets. For < k < d we denote by fk(V) the number of fc-faces of 
V . Note that every fc-face F of V is also a /c-polytope whose faces are all the faces of V contained 
in F . A fc-simplex in K d , k < d, is the convex hull of any k + 1 affinely independent points in K d . 
A polytope is called simplicial if all its proper faces are simplices. Equivalently, V is simplicial if 
for every vertex v of V and every face F € V , v does not belong to the affine hull of the vertices in 
F\{v}. 

A polytopal complex C is a finite collection of polytopes in K d such that (i) £ C, (ii) if V £ C 
then all the faces of V are also in C and (iii) the intersection V n Q for two polytopes in C is a face of 
both V and Q. The dimension dim(C) of C is the largest dimension of a polytope in C. A polytopal 
complex is called pure if all its maximal (with respect to inclusion) faces have the same dimension. 
In this case the maximal faces are called the facets of C. We will use the term d-complex to refer to 
a pure polytopal complex whose facets are ci-dimensional. A polytopal complex is simplicial if all its 
faces are simplices. Finally, a polytopal complex C is called a subcomplex of a polytopal complex 
C if all faces of C are also faces of C. 

One important class of polytopal complexes arise from polytopes. More precisely, a (i-polytope 
V , together with all its faces and the empty set, form a polytopal cf-complex, denoted by C(V). 
The only maximal face of C(V), which is clearly the only facet of C(V), is the polytope V itself. 
Moreover, all proper faces of V form a pure polytopal complex, called the boundary complex C{dV). 
The facets of C{dV) are just the facets of V , and its dimension is dim('P) — 1 = d — 1. Given a 
polytope V and a vertex v of V , the star of v is the polytopal complex of all faces of V that contain 
v, and their faces. The link of v is the subcomplex of the star of v consisting of all the faces of the 
star of v that do not contain v. 

The /-vector (/_i('P), fo(V), . . . , fd-i(V)) of a simplicial (i-polytope V is defined as the (d+ 1)- 
dimensional vector consisting of the number fk(V) of A;-faces of V , —1 < k < d, where f_i(V) = 1 
refers to the empty set. The /i-vector (ho(P), hi(V), ■ ■ . , hd(V)) of a simplicial ci-polytope V is 
defined as the (d + l)-dimensional vector, where h k {V) := Ei=o( _1 ) fe_ * {t^fi-X^P), < k < d. 
The number hk(V) counts the number of facets of V in a shelling of V , whose restriction has size 
k; this number is independent of the particular shelling chosen (cf. |28[ Theorem 8.19]). It is easy 
to verify from the defining equations of the ^('P)'s that the elements of the /-vector determine the 
elements of the ii-vector and vice versa. Moreover, the elements of the /-vector (or, equivalently, the 
h- vector) are not linearly independent; they satisfy the so called Dehn-Sommerville equations, which 
can be written in a very concise form as: h^V) = hd-kiP), < k < d. An important implication 
of the existence of the Dehn-Sommerville equations is that if we know the face numbers fk(V) for 
all < k < |_f J — 1, we can determine the remaining face numbers fk(P) for all |_f J < k < d — 1. 

In what follows we recall some facts from |17[ Section 5.2] that will be of use to us later. Let V 
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be a d-polytope in K d , F a facet of V, and H the supporting hyperplane of F (with respect to V). 
For an arbitrary point p in K d , we say that p is beyond (resp., beneath) the facet F of V, if p lies in 
the open halfspace of H that does not contain V (resp., contains the interior of V). Furthermore, 
we say that an arbitrary point v' is beyond the vertex v of V if for every facet F of V containing 
v, v' is beyond F, while for every facet F of V not containing v, v' is beneath F. The vertices of 
the polytope V' = CH&({P\{v~Y) U {«'}) are the same with those of V, except for v which has been 
replaced by v' . In this case we say that V' is obtained from V by pulling v to v'. The vertex v' of 
V' does not belong to the affine hull of the vertices in F'\{v'} for every face F' of V 1 . The following 
result is well-known. 

Theorem 1 ([20}l22]). Let V be a d-polytope. 

(i) The d-polytope V' we obtain by pulling a vertex ofV has the same number of vertices with V , 

and f k (T) < f k (P') for all 1 < k < d - 1. 
(ii) The d-polytope V' we obtain by successively pulling each of the vertices ofV is simplicial, has 
the same number of vertices with V , and fk(P) < fk(T~") for all 1 < k < d — 1. 

In the rest of the paper, when we refer to parallel hyperplanes we assume that they have the 
same unit normal vector, i.e., they have the same orientation. Moreover, if two hyperplanes II and 
II' are parallel, we say that II' is above IT if II' lies in the positive open halfspace delimited by IT. 

Let & = {Vi,T > 2, • • • , Vm\ be a set of m ci-polytopes lying on m parallel hyperplanes IIi, H2, . . . , 
n m of E d+1 , respectively. Throughout this section we assume that m > 2 is fixed, and that IF, is 
above IT for all j > i. We denote by Pj the set of vertices of Vi, by 7ij the cardinality of Pi, and 
by P the union P = P\ U P2 U . . . U P m . Let V = CH ( i + i(P); note that, for each i, not all vertices 
in Pi are necessarily vertices of V . Furthermore, among the polytopes in only V\ and V m are 
faces of V . 

The theorem that follows is the adaptation of Theorem [T] in the context of m parallel polytopes. 
Again, we want to perturb the points in P so that V' (the polytope we obtain after perturbing 
the points in P) is simplicial with fk(V) < fk(P') f° r k > 1, but we want to retain the property 
that the points lying on a hyperplane ITj, if perturbed, are replaced by points that lie on the same 
hyperplane. This is almost possible. More precisely, all the faces of the polytope V 1 are simplicial, 
with the possible exception of V[ = V' Pi IIi and V' m = V D H m . 

Lemma 2. Let & = {V\,V2, ■ ■ ■ ,P m } be a set ofm > 2 d-polytopes lying on m parallel hyperplanes 
IIi,ri2, . . . ,II m o/E rf+1 , respectively, where ITj is above IT for all j > i. Let Pi be the vertex set of 
Pi, 1 < i < m, P = P\ U P2 ■ ■ ■ U P m , and V = CHd+i{P)- The points in P can be perturbed in 
such a way that: 

(i) the points of P in each hyperplane IT remain in Hi, 1 < i < m, 

(ii) all the faces ofV', except possibly the facets V[ and V' m , are simplices, and, 
(in) f k {V) < f k {V) for all 1 < k < d, 

where V' is the polytope we obtain after having perturbed the vertices ofV in P. 

Proof We construct V' in three steps. Firstly, we properly perturb the points in IIi so that V[ is 
simplicial, then we do the same for the points in II m , and, finally, we pull every vertex of V in II; 
for 2 < i < m — 1. 

Let v G IIi and choose any v' in IIi beyond all the facets of V other than V\ that contain v. 
If we focus our attention on the hyperplane ITi and on the polytope V\ G ITi we have that v is a 
vertex of V\ and v' is some point beyond all the facets oiV\ containing v. We consider the polytope 
V' = CH d+1 ((P\{v}) U {v'}). In view of pH Theorems 5.1.1 & 5.1.2] the £>faces of V' are either: 
(i) A;- faces of V[, or 
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(ii) fc-faces of P not in III that do not contain v, or 

(iii) faces of the form CHd+i({v'} U Gk-i), where Gk-i is a (/c — l)-face not containing v of a facet 
P III of P containing v. 

The number of fe-faces of V' does not change due to Case (ii). The number of fe-faces of V' may 
increase in Case (iii) if G^-i is contained in a non-simplicial fc-face of P, otherwise it does not 
increase. Finally, Case (i) is a consequence of Theorem [T]^i) . 

We set P := V and we repeat the above procedure for every point v £ P\. After having 
perturbed all points in Pi we obtain a polytope P' such that P( is a simplicial polytope. Moreover, 
we have the additional properties that /fc(P) < /fc(P')j f° r k > 1, and that for every v' £ Iii and 
for each face F 1 £ V' not lying in 111, the point v' does not lie in the affine hull of the vertices in 
F'\{v'}. We repeat exactly the same procedure on the points of H m . Again, we denote by P the 
polytope we obtain after having perturbed the points in Pi and P m . 

It is left to pull the points in the remaining hyperplanes; in fact, we need only perturb the points 
in P that are vertices of V . Consider a vertex v of P, such that v £ Pi for some 2 < i < m — 1. 
Let v' be any point in IIj that lies beyond all the facets of V containing v. The choice of such a 
v' is possible since Ilj is a hyperplane containing v but not a supporting hyperplane of V . The 
polytope V' = CHd+i((P\{v }) U {v'}) is the one we obtain by pulling the vertex v to v' and thus 
from Theorem [T^i) we have /fc(P) < /fc(P') for all k > 1. We continue the same procedure for 
every vertex v of V in Ilj, and 2 < i < m — 1. After having pulled all the vertices of V in Ilj, for 
2 < i < m — 1, we obtain a polytope P' with the property that for every vertex v' of P' in Ilj, with 
2 < i < m — 1, and for every face P' £ P', the point i/ does not lie in the affine hull of the vertices 
in F'\{v'}. 

Summarizing all of the above, we deduce that, after having perturbed all the vertices of P in 
P, we get a polytope P' with the same number of vertices as P, such that: (1) /fe(P) < /fc(P') for 
all k > 1, (2) PJ and P^ are simplicial, and (3) for every vertex v' of P' and every face P' G P', 
such that P' does not lie in LTi or II m , the vertex v' does not lie in the affine hull of the vertices in 
F'\{v'}. The latter implies that all the faces of P', not in 111 or Il m , are simplices. This completes 
our proof. □ 

In view of Lemma [2j it suffices to restrict our attention to sets of polytopes where P is 
simplicial with the possible exception of its two facets Pi and P m . Let II be any hyperplane 
between and parallel to the hyperplanes II m _i and II m and consider the intersection P := P n II 
(see Fig. [T]). Let J- be the set of faces of P having non-empty intersection with II. Note that P is a 
d-polytope, which is, in general, non-simplicial, and whose proper non-trivial faces are intersections 
of the form P Pi II where F G P. Let A = (ai, a2, • • • , a m ) and B = (f3±, fa, ■ ■ ■ , /3 m ) be two vectors 
in N m . We say that A =4 B if aj < /3j for all 1 < i < m, and denote by \A\ the sum of the elements 
of A, i.e., |^4| = a i- The following lemma provides an upper bound on the number of fc-faces 

of P. 

Lemma 3. The number of k-faces of P is bounded from above as follows: 



in 



(0,...,0,l)^A^!(fe,...,fe) i=l 
\A\=k+\ 

where aj is the i-th coordinate of the vector A £ N m . 

Proof. According to Lemma [2j it suffices to consider the case where P is simplicial except possibly 
for its facets Pi and P m . In this context, a fc-face P £ P is simplicial and it is defined by k + 1 
vertices of P, where at least one vertex comes from P m , whereas the remaining k vertices are vertices 
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Figure 1: An example of four polytopes lying on four parallel hyperplanes of E . The white vertices 
are not vertices of V . The polytopes V2 and V3 are shown in black. All faces of C are shown in 
green (only the facet V\ is shown), whereas all faces of V± are shown in red. The faces in blue are 
faces of P, whereas the faces in brown are faces of V = V Pi II. 

of P±,P2, . . . , Pm-\. Let cti be the number of vertices of F from Pi. Clearly, we have < oti < k 
for 1 < i < m — 1, 1 < a m < k, and a\ + ai + . . . + a m = k + 1. The maximum number of possible 
(«i — l)-faces of Vi is ( „ ), which implies that the maximum possible number of /c-faces of T is 

lli=i Summing over all possible values for the ctj's we get the desired expression. □ 

Using the bounds for the face numbers of T from Lemma [3l we arrive at the following lemma 
concerning bounds on the number of /c-faces of V, for small values of k. 

Lemma 4. Let ni = fo(Vi), 1 < i < m. The following asymptotic bounds hold: 

m 

fk (v)=o( £ +£-r m{fc+1 ' bJ} ), o<k<n\. 

l<ij^j<m i=l 

Proof. In what follows k will be at most |_f J , while we denote by k! the quantity k' = min{£;+l, J }■ 
Let T{m) be the worst-case complexity of fk(V). We are going to prove, by induction on m, that 
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for all m > 1, T(m) < c{^2, l<i ^L- <m nin^ + n i)> where c is some appropriately large constant 

that depends only on d. 

The case m = 1 is trivial since the number of fc-faces of a (i-polytope with n\ vertices is in 
0(n\ ). Let us now assume that m > 2 and that our statement holds for m — 1; we shall prove it 
for m. To this end, we consider a set of m <i-polytopes & = {Vi,T , 2, ■ ■ • ,P m }, lying on m parallel 
hyperplanes III, . . . , Tl m of E rf+1 , such that Tlj is above Ilj for all j > i. Let m = foi'Pi), 1 < i < m, 
and denote by V the convex hull CH<i+i{&)- Consider a hyperplane IT parallel to Ti m -\ and Yi m 
and between them, and let V = V fl II (refer to Fig. [I]). We denote by £ the set of faces F of V 
such that F P m and F n ft = 0, and by J 7 the set of faces F oiV with F n ft / 0. Clearly, the 
set of faces of V is equal to the disjoint union of C, J- and the set of faces of V m \ hence: 

fk(V) = f k (C) + f k (J) + HV m ). (1) 

By the induction hypothesis we have that the number of fc-faces of CH c i + i(£P\{'P m }) is at most 
T{m — 1). Since the fc-faces in C are fc-faces of C H ( i+i(r0 g \{'P m }) , we have that fk(£) is at most 
T(m — 1). On the other hand, by Lemma [3j we immediately get that fk{J~) = 0{n^ n Y^ 7 ^ 1 n % + 
rimY^i^i 1 n i)- Finally, since V m is a d-polytope, fk(P m ) = O(nJ^). Combining these bounds with 
eq. dl]), we arrive at the following recurrence relation for T{m): 

m—X m—1 

T{m) < T{m - 1) + 0{n k m ^ m + n m n i) + °( n ™)- 

i=l i=l 

It is straightforward to verify that T(m) satisfies: T(m) < c (X/i<»y<m n i n j + S^=i n i )> ^ or some 
appropriately large constant c (that depends only on d); this establishes our claim. □ 

Exploiting the bounds from Lemma H] for < k < Nfj, and using the Dehn-Sommerville equa- 
tions of an appropriately defined simplicial (d+ l)-polytope containing all faces in V, except for the 
facets V\ and V m , we derive asymptotic bounds on fk{T > ), for all L^y^J < k < d. Our results are 
summarized in the following lemma. 

Lemma 5. Let rii = fo(Vi), 1 < i < m. The following asymptotic bounds hold: 

fk(V)=0( ]T n t nf V ), [d±k\<k<d. 

Proof. Let y (resp., z) be a point below LTi (resp., above II m ), such that the vertices of V\ (resp., 
Vm) are the only vertices of V visible from y (resp., z). To achieve this, we choose y (resp., z) to be 
a point beyond the facet V\ (resp., V m ) of V, and beneath every other facet of V. Let Q be the set 
of points consisting of y, z and the vertices of V, and let Q = CHd+i{Q) (refer to Fig. [2]). Observe 
that the faces of V, except for the facets V\ and P m , are all faces of Q. To see that, notice that a 
supporting hyperplane Hp for a facet F G V, with F ^ Vi,V m , is also a supporting hyperplane for 
Q. Indeed, the vertices of F are vertices in Q \ {y, z} = P, and thus, every vertex in P that is not 
a vertex of F strictly satisfies all hyperplane inequalities for V . Also, by construction, the points y 
and z strictly satisfy all hyperplane inequalities apart from those for E4 and II m , respectively. Since 
Hp is a hyperplane other than III and II m we deduce that all points in P U {y, z} lie on the same 
halfspace defined by Hp, and therefore Hp supports Q. The faces of Q that are not faces of V are 
the faces in the star S y of y and the star S z of z. To verify this, consider a fc-face F of "Pi, and let 
Fx be a face in V\ that contains F. Let H\ be a supporting hyperplane of F\ with respect to V. 
Tilt Hi until it hits the point y, while keeping H\ incident to F\, and call H[ this tilted hyperplane. 
H[ is a supporting hyperplane for y and the vertex set of V±, and thus is a supporting hyperplane 
for Q. The same argument can be applied for the star of z. In fact, the boundary complex dV\ 
(resp., dV m ) of V\ (resp., V m ) is nothing but the link of y (resp., z) in Q. 



S 



z 




It is easy to realize that, for k < d, the set of k- faces of Q is the union of the /c-faces of V, S y 
and S z , where we have double-counted the /c-faces of dV\ and &P m . For k = d the set of facets of 
Q is the union of the set of facets of S y , S z and V, except for its two facets V\ and V m . We can 
concisely write the relations described above as: 

fk(Q) = fk{V) + f k {Sy) + fk{S z ) ~ fk{Vi) - fk(Pm), 0<k<d. (2) 

where fd(Vi) = fdi'Pm) = 1- The /c-faces of Q in S z are either /c-faces of &P m or /c-faces defined by 
z and a (k — l)-face of dV m - In fact, there exists a bijection between the [k — l)-faces of &P m and 
the /c-faces of S z containing z. Hence, we have: 

fk(S z ) = f k (&P m ) + fk-i(8P m ), 0<k<d, (3) 

where f-i(&P m ) = 1 and fd,(dV m ) = 0. Analogously, the fc-faces of Q in S y are either fe-faces of 
dV\ or /c-faces defined by y and a (k — l)-face of dV\. As for S z , there exists a bijection between 
the (k — l)-faces of dV\ and the /c-faces of S y containing y. Hence, we have: 

fk(S y ) = f k (m) + fk-iidVi), 0<k<d, (4) 

where f-^dVi) = 1 and fd(&P\) = 0. Therefore, relation ([2]) can be rewritten as: 

fk(Q) = fk{V) + f k (8Pi) + h-iim) - fk{Vx) + fk(9Pm) + fk-i{dV m ) - f k (V m ), < k < d. 

or, in the much simpler form: 

fk{Q) = fk(V) + /fc-i(SPi) + fk-i{dV m ) - a k , 0<k<d, (5) 

where a k = 0, for < k < d, whereas = 2. 

Let us now turn our attention on deriving bounds for the face numbers fk(Q)- Since dVg, 
£ = l,m, is the boundary complex of the n^-vertex ci-polytope Ve, we have: 

hi&Pt) = 0(nf n{k+1 '^ 1} ), 0<k<d, e = l,m. (6) 

Combining the bounds from Lemma [4] with relations ([5]) and ([6]) , we get the following bounds for 
the lower-dimensional face numbers of Q: 

m 

fk(Q) = 0{ £ ^ + £nf n{fc+1 ' L5J} ) = 0( 0<fc<L|j. (7) 

l<i7^j<m i=l l<*7^i<TJ 

However, Q is a simplicial {d+ l)-polytope, since the facets of V are simplicial, and since all facets 
of S y and S z are simplicial (the facets of S y and S z are defined via simplicial id — l)-faces of V 
and the points y and z, respectively). Let us recall the defining equations for the elements of the 
h- vector of Q in terms of the elements of the /-vector of Q: 

k 

hk{Q) = Y(-i) k - l ( d l\~_ l ^)h-i{Q), 0<fc<d + l. (8) 

i=0 

Combining equations ([8]) with relations ([7]), as well as the fact that /_i(Q) = 1, we get: 

h k (Q)=0( Yl n *4 lJ )' 0<k<[^\. (9) 

l<i^j'<m 

We are now going to use the Dehn-Sommerville equations for Q to bound the number of /c-faces 
fk(Q) °f Qi f° r k > L^T^J- The Dehn-Sommerville equations can be rewritten as follows (cf. |28[ 
Section 8.4]): 

d+l 

A-<« = E'(("^70 + G-/i +i ))wa>- l^j<*<*+i, a") 
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6 
2 

denotes the sum where the last term is halted if and only if the quantity | 

i=0 

is integral (which is our case since in our setting 5 = d+1, which is even). Combining relations ([9|) 
and (|10p . we get: 

fk-i(Q) = o( Yl ^ n | fJ )> L^J<*<d + i. (ii) 

l<i^j'<m 

Now using relation ((5|), we arrive at the following bounds on the number of A:-faces of V for \ ^-^-\ < 
k < d: 

fk(V) = f k (Q) - f k -i{&P x ) - f k -i(dV m ) + a k < f k {Q) = C( E n * n i fJ )' 

where we used the fact that fd-iidVi) > ad, since fd-iidVi) > d+ 1 > 4. □ 
By Lemmas 2] and O we deduce that the worst-case complexity of the convex hull CHd+\(&) 

is 0(^i<j^j<m n i n j 2 )• A s we wm see m Subsection 14.21 (see Corollary 1 1 1 [) . this bound is asymp- 
totically tight for any odd d > 3. Hence: 

Theorem 6. Let & = {Vi,T , 2, ■ ■ ■ , "P-m} be a set of a fixed number of m > 2 d-polytopes, lying on 
m parallel hyperplanes o/E d+1 ; where d>3 and d is odd. The worst-case complexity of the convex 

hull CHd+ii^) is 6(Ei<,^< m n t"j 2 )> where m = fo{Vi), I < i < m. 

3 Computing the convex hull of parallel convex polytopes 

The upper bound on the worst-case complexity of CHd+i(&) in Theorem suggests that, in order to 
compute CHd+\{&) for d > 3 odd, it pays off to use an output-sensitive algorithm for constructing 
the convex hull of the point set P formed by the vertices of the TVs. 

Let us briefly discuss the output-sensitive algorithms applicable in our setting; in what follows 
ni is the number of vertices of Vi, n is the total number of input points and / the size of the output. 
The dimension d + 1 > 4 below is the dimension of the ambient space E rf+1 and is considered to 
be even; the polytopes Vi are then <i-dimensional. One of the earliest algorithms is Seidel's shelling 
algorithm |27| that runs in 0(n 2 + /logn) time. The preprocessing step of Seidel's algorithm 
was later on improved by Matousek and Schwarzkopf [21] . resulting in an 0(n 2_2 ^L( rf + 1 )/ 2 J+ 1 )+ e _|_ 
/logn) time algorithm, for any fixed e > 0. Chan, Snoeyink and Yap describe a divide- 
and-conquer algorithm for constructing four-dimensional convex hulls in 0((n + /)log 2 /) 

time. 

Finally, Chan |10| improved the gift- wrapping algorithm of Chand and Kapur |12| . yielding an 
0(n log / + (n/) 1 -V(L(d+i)/2J+i) log o(i) ^ time a i gorit hm. 

In our case we have / = 0(^ 1<i ^- <m n^n" -1 ), where a = Y^rl > 2 (recall that we consider 
only odd d). Applying the algorithms above we get the following. Seidel's algorithm runs in 
0(n 2 + (X/i<i^j<m n i n< j~ l ) 1°§ n ) ti me - Matousek and Schwarzkopf's improvement results in an 
0(n 2_2/ '( a+1 * >+<E + (X^i<i^<m n i n< j~ 1 ) l°g n ) time algorithm. Chan's algorithm yields a running 
time of 0((X^i<j^j<m n r^ 0+1 ^ n J ^ a+1 ^) log *- 1 ) n). Finally, the algorithm by Chan, Snoeyink and 
Yap gives 0((^i<j^j<m n i n j) l°g 2 n )- F° r an d > 3, i.e., for a > 2, Chan's algorithm is never faster 
than the other two algorithms. For d > 5, i.e., a > 3, it makes no difference to choose between the 
algorithm of Seidel and that of Matousek and Schwarzkopf. For d = 3 the situation is a bit more 
complicated. For some choices for the n^'s Matousek and Schwarzkopf's algorithm outperforms all 
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other algorithms, whereas for other choices of the rij's the algorithm of Chan, Snoeyink and Yap is 
the best choice. Consider, for example, the cases: 

(1) rt\ = 0(n), ?i2 = ) an( i n i = 0(1) for 3 < i < m, and 

(2) ni = 6(n), n 2 = and n* = 0(1) for 3 < t < m. 

In the first case Chan, Snoeyink and Yap's algorithm runs in o(ra 4 / 3+e ) time and outperforms all 
other algorithms. In the second case, Matousek and Schwarzkopf's algorithm seems to be the best 
choice, since its running time is 0{ri2n\ogn) (compared to Chan, Snoeyink and Yap's algorithm, 
the running time of which is 0(n2nlog 2 n)). Summarizing, we arrive at the following: 

Theorem 7. Let = {Vi,T > 2, ■ ■ ■ ,V m } be a set of m d-polytopes, lying on m parallel hyperplanes 
ofK d+1 , where d > 3 and d odd. Let m = fo(Vi), 1 < i < m, and n = Y^iLi n i- We can compute 

the convex hull CH ( i+i(& > ) in 0((Xa<i=y<m n i n j 2 )l°g n ) time for d > 5, and 0(min{n 4//3+e + 
(X)l<t#<m n i n i) lo § n > (Ei<i^< m n i n i) lo S 2 n i) time f or d =3> f or an V fi xed e > 0. 

4 Convex hulls of spheres with a fixed number of distinct radii 

In this section we derive tight upper and lower bounds on the worst-case complexity of the convex 
hull of a set of spheres in E^ having m distinct radii, where m is considered to be fixed. 

Let S be a set of n spheres o~k = (cf., r^), 1 < k < n, in E d , and let CiJrf(S) be the convex hull of 
the spheres in E. A face of circularity £ of CiJd(S), < £ < d — 1, is a maximal connected portion 
of the boundary of Cffd(E) consisting of points where the supporting hyperplanes are tangent to a 
given set of (d — £) spheres of S. In the special case where all spheres have the same radius, CHd(T,) 
is combinatorially equivalent to the convex hull CHd(K) of the centers K of spheres in S, in the 
sense that each face of circularity £ of CH^i^) corresponds to a unique (d — i— l)-face of CHd(K), 
for < I < d - 1. 

We consider here the case where the radii r^ can take m distinct values, i.e., £ {pi,P2, ■ ■ ■ > Pm}- 
Without loss of generality we may assume that < p\ < p 2 < ■ ■ ■ < p m - We identify E d with the 
hyperplane Hq = = 0} of E rf+1 and we call the {d + l)-axis of E d+1 the vertical axis, while 

the expression above will refer to the (d + Incoordinate. Let IL;, 1 < i < m, be the hyperplane 
= pi} in E rf+1 , and let P be the point set in E d+1 obtained by mapping each sphere to the 
point pk = (cfc,rfc). Let Pi denote the subset of P containing points that belong to the hyperplane 
Ilj, and let rij be the cardinality of Pi. We denote by V the convex hull of the points in P (i.e., 

V = CHd+i(P))- We further denote by Vi the convex hull of the points in Pi (i.e., Vi = CHd(Pi)); 
more precisely, we identify Ilj with E rf , and then define Vi to be the convex hull of the points in 
Pi, seen as points in E d . We use to denote the set of the ViS. Let Pi be the subset of Pi that 
defines Vi (i.e., the points in Pj are the vertices of Vi and thus Vi = CH^Pi)), and let fii < Hi be 
the cardinality of Pi. Finally, let P = (J2=i Pi an d V = CH ( i + i(P). Notice that it is possible that 

V 7^ V, such a situation will arise if Pi ^ P\ (resp., P m ^ P m ), in which case the intersection of V 
with ITi (resp., II m ) will consist of more than one <i-face of V. On the other hand V and V have 
the same interior. 

Let Ao be the half lower hypercone in E d+1 with arbitrary apex, vertical axis, and angle at the 
apex equal to j. By X(p) we denote the translated copy of Ao with apex at p; observe that the 
intersection of the hypercone \(pk) with the hyperplane Hq is identical to the sphere cifc. Let A be 
the set of the lower half hypercones {A(pi), X(p2), ■ ■ ■ , A(p n )} in E d+1 associated with the spheres of 
S. The intersection of the convex hull CHd+i(A) with Hq is equal to CHd(T<). 

Let O' be a point in the interior of V. We then have the following: 
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Theorem 8 ([61 Theorem 1]). Any hyperplane of E rf supporting CiJd(E) is the intersection with 
Hq of a unique hyperplane H of E d+1 satisfying the following three properties: 

1. H supports V , 

2. H is the translated copy of a hyperplane tangent to Xq along one of its generatrices, 

3. H is above O' . 

Conversely, let H be a hyperplane ofE d+1 satisfying the above three properties. Its intersection with 
Ho is a hyperplane o/E d supporting Cif^(S). 

Theorem [8] implies an injection ip : CHd(T,) — > V that maps each face of circularity (d — I — 1) 
of Ci?rf(E) to a unique £-face of V, for < £ < d— 1. Theorem [8] also implies that points in Pi\Pi, 
1 < i < m, can never be points on a supporting hyperplane H of V satisfying the three properties 
of the theorem. Therefore, if is, in fact, an injection that maps each face of circularity (d — I — 1) of 
Cifd(E) to a unique £-face of V , < I < d — 1. Observe that V is the convex hull of the set & of 
m convex polytopes lying on m parallel hyperplanes of E d+1 . By employing Theorem [6] of Section 

\d\ \—\ 

O we deduce that "P's complexity is 0{^ x<i ^ <rn n i nj' z ) = 0(J2i<i^j< m n i n j 2 )> which, via the 
injection tp : CHd(T<) — > V, is also an upper bound for the worst-case complexity of CHd(T,). 

4.1 Lower bound construction with two distinct radii 

For any even dimension 26, the trigonometric moment curve 72^(t) in E 2<5 is the curve: 

72^ (t) = (cost, sin t, cos 2i, sin 2i, ... , cos <5i, sin t£ [0, tt). 

Notice that points on 7^(t) are points on the sphere of E 2<5 centered at the origin with radius equal 
to \f5. For any set P of n points on 7^ (t), the convex hull CH 2 s(P) is a polytope Q combinatorially 
equivalent to the cyclic polytope C2s(n) (cf. [13 [28]). Therefore, /2<5-i(Q) = &(n s ). 

Suppose now that the ambient space is ¥, d , where d > 3 is odd. Let Hi be H 2 be the hyperplanes 
{xd = zi} and {xj = z 2 }, where zi,z 2 G K and z 2 > z\ + 2(?i2 + 2)\^5; the quantity n 2 will be 
defined below. Consider a set Si of n\ + 1 points, treated as spheres of E d of zero radius, on 
the [d — l)-dimensional trigonometric moment curve 7^!li(t) embedded in (please refer to Fig. 
|3](left) , as well as Fig. |4] for the view of the construction from the positive a^-axis). Among the 
n\ + 1 points, the first n\ points are chosen with t G (0, ^), whereas for the remaining point we 
require that t S ,7r). This implies that the xi-coordinate of the first n\ points of Si is positive, 
whereas the xi-coordinate of the last point of £1 is negative. Let Y> 2 be the projection, along the 
x^-axis, of Ei on the hyperplane H 2 . Clearly, the n\ + 1 points of T, 2 in H 2 lie on the (d — 1)- 
dimensional trigonometric moment curve r y^ r _ 1 (t) embedded in H 2 . The points of Ej, i = 1,2, lie on 
a (d — 2)-dimensional sphere of E d , centered at the point (0, 0, ... , 0, Zi), with radius \fb. Moreover, 
the number of facets of the polytope Qj = CH^_i(T,i) is Q(n 1 2 ) = G(n 1 2 ). The convex hull 

of the 2(ni + 1) points of Ei U £2 is a prism A. A consists of 0(n 1 2 ) facets not lying on H\ or 
H 2 , called the vertical facets of the prism, the (d — l)-face of CH^iiTii), called the bottom facet, 
and the (d — l)-face of CHd_i(E 2 ), called the top facet. For each vertical facet F of A, we denote 
by vp the unit normal vector of F pointing outside A, and by F + (resp., F~) the positive (resp., 
negative) open halfspace delimited by the supporting hyperplane of F. Regarding the ridges of A, 
those that are intersections of vertical facets of A will be referred to as vertical ridges. Notice that 
the vertical ridges of A are perpendicular to H\ and H 2 . 

Let Y be the oriented hyperplane {37 = 0} with unit normal vector V = (1,0, ... ,0). Let also 
Y + and Y~ be the positive and negative open halfspaces of E d delimited by Y, respectively. Y 
contains the a^-axis, and is perpendicular to the hyperplanes Hi and H 2 . Recall that n\ points 
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of Ej are contained in Y + , whereas exactly one point of Ej is contained in Y~. Clearly, Y is in 
general position with respect to A, Q\ and 0,2- Let Qi be the intersection of with Y, and let 
Ti be the set of faces of Qi intersected by Y. Qi is a (d — 2)-polytope, and its number of vertices 
is at most m, since Y cuts at most ri\ edges of Qi. This implies that the complexity of Qi is 

i d ~ 2 j I -J— i 

0(n 1 2 ) = 0(n 1 2 ); the same bound holds for Ti. Since there are no facets of A in Y~ (Y 

L-J 

contains a single vertex of Qi), and since the number of facets of Qi is 0(n 1 2 ), we conclude that 

L-J 

the number of facets of Qi contained in Y + is also 0(n 1 2 ); the same bound holds for the number 
of vertical facets of A in Y + . 

Define now a set E3 = {ao,ai, . . . ,cr n2 +i} of ri2 + 2 spheres in E d , where Ok = (c^,/?), and 
c k = (0, . . . , 0, (2k + 1)^/6), < k < ri2 + l. In other words, the sphere <7fc is centered on the x^-axis, 
with the d-th. coordinate of its center Ck being (2k + l) - ^, while its radius is p. We choose p to be 
smaller than \J~S, but large enough so that each sphere Oi satisfies the following two properties: 

(1) it does not intersect any of the ridges of A (including the vertical ridges of A), and 

(2) it intersects the interior of all vertical facets of A. 

Notice also that for this choice for p, none of the spheres in E3 intersects the hyperplanes H\ and 
H2 (recall that 22 > z i + 2(ii2 + 2)y/S), while the spheres in E3 are pairwise disjoint; these two 
observations, however, are not critical for our construction. 

We are now going to perturb the centers of the spheres in E3 to get a new set of spheres E 3 (see 
Fig. fright), as well as Fig. [5] for the view of the construction from the positive x^-axis). Define a' k 
to be the sphere with radius p and center c' k = + (X^ =0 2^)^ = c k+ e (2 — ^r)^ 5 where 0<£<1. 
The quantity e is chosen so that the spheres in E 3 satisfy almost the same conditions as the spheres 
in E3. In particular, we require that condition (1) is still satisfied, while we relax the requirement 
on condition (2): we now require that a' k intersects the interior of all vertical facets of A contained 
in Y + . In addition to the two conditions above, we also require that for each k, < k < n-i + 1, the 
(d — 2)-dimensional sphere Ok n a' k is contained in F~ for all vertical facets F of A in Y + . 

We will now show that for each pair (a' k , F) , where 1 < k < rii and F is a vertical facet of A in 
Y + , the spherical cap F + Pi a' k induces a face of circularity [d — 1) in Ci?d(E). Let F\ and F2 be the 
ridges of A on the boundary of F contained in the top and bottom facet, respectively. Finally, let 
Sk be the supporting hyperplane of parallel to F; we consider Sk to be oriented as F (i.e., the 
unit normal vector of Sk is z/p), and thus Ok lies in the closure of the negative halfspace delimited 
by Sk- Notice that Sk is also a supporting hyperplane for E3. Let S' k be the hyperplane we get by 
translating Sk by the vector e(2 — -^)v. S' k supports a' k , but fails to be a supporting hyperplane 
for E 3 . More precisely, S' k intersects all spheres a'j with j > k, whereas all spheres a'j with j < k, 
are contained in the negative open halfspace delimited by S' k . We can, however, perturb S' k so that 
it supports E 3 : simply slide S' k on sphere a' k towards F±, while maintaining the property that it 
remains parallel to F\ and Fi- We keep sliding S' k until it has empty intersection with any sphere 
u'a with j > k. Notice that due to the way we have perturbed the centers of the spheres in E3 to 
get S 3 , the new hyperplane S'l we get via this transformation is a supporting hyperplane for E 3 . 
In fact, S'l is a supporting hyperplane for the sphere set E = Ei U E2 U E 3 (it touches C//rf(E) at 
a' k only), which implies that S k corresponds to a unique face of circularity (d — 1) on Cif^(E). 

The same construction can be done for all k with 1 < k < ri2, and for all vertical facets of 

L~J L~J 
A in Y + . Since we have B(n x 2 ) vertical facets of A in Y + , we can construct n2©(^ 1 2 ) distinct 

supporting hyperplanes of CiJrf(E), corresponding to distinct faces of circularity (d— 1) on CH^iTi). 

L-J 

Hence the complexity of CH^(Tt) is J7(n2n 1 2 ). Without loss of generality, we may assume that 

\ — \ [— J I — J 

fi2 < ni, in which case we have n.2?^ 2 > ^(ri2?T- 1 2 + nin 2 2 )• Hence, we arrive at the following: 
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Figure 3: The lower bound construction in the case of two radii. The points in Si (resp., S 2 ) are 
shown in black (resp., white). The hyperplane Y is shown in green, while the prism A is shown in 
black. The facet F in blue is one of the vertical facets of A in Y + . The sphere sets S3 (left) and 
S 3 (right) are shown in red. The red spherical caps on the left correspond to a unique supporting 
hyperplane of CHd{^i U S2 U S3). The red spherical caps on the right correspond to faces of 
CiJ d (Si U S 2 U S' 3 ) of circularity (d - 1). 
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Figure 5: View from the positive x^-axis of the construction in Fig. fright). 
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Theorem 9. Fix some odd d > 3. There exists a set £ of spheres in E , consisting of n% spheres of 

L-J 

radius p^, i = 1,2, u/ii/i pi 7^ p2» suc/i i/iai i/ie complexity of the convex hull CHd(T>) is Q(nin 2 2 + 

LfU 
n 2 re 1 / J. 

4.2 Lower bound construction with m distinct radii 

We can easily generalize the lower bound construction of the previous subsection in the case where 
we have rii spheres of radius pi, 1 < i < m, m > 3, and the radii pi are considered to be mutually 
distinct. 

As in the previous subsection, the ambient space is E rf , where d > 3 is odd. Let N\ = n « 
and iV 2 = n\. We construct the set S = Si U S2 U S 3 as in the previous subsection where Si and 
S2 contain each N\ + \ points and S 3 contains iV 2 + 2 spheres of some appropriate radius p (recall 

that in the construction of the previous subsection p ~ \J dj ^- > 1). We then replace rii among the 
N± points of Si (resp., S2) contained in Y + by spheres with the same center and radius equal to 
r l , where < r < 1. Furthermore, we replace the unique point of Si (resp., £2) hi Y~ by a sphere 
of the same center and radius r 2 . We choose r small enough so that the following two conditions 
hold: 

(1) the prism A r = Cff^Si US2) is combinatorially equivalent!! to the prism Ao (this is the prism 
we get for r = 0, which is the prism A of the previous subsection), and 

(2) the two requirements for the spheres in S 3 are still satisfied: each a' k does not intersect any of 
the ridges of A r , while each a' k intersects the interior of all vertical facets of A r contained in 
Y+. 

As described in the previous subsection, the convex hull CHd(Ti) of the set S = Si U S2 U S 3 of 

L-J 

2(iVi + 1) + N2 + 2 spheres has A r 2G(A^ 1 2 ) faces of circularity (d — 1), and hence its complexity 
is Jl(iV2-/v}^ ) = ^(rai(^™ 2 n i)^^)- Without loss of generality we may assume that ri2 > n± > rij 
for all 3 < i < m, in which case we have: ni(^™ 2 m) ^J > mn 2 2 > m ^_^ (52l<i^j<m n i n j 2 )■ 

L-J 

Since m is fixed, we conclude that the complexity of CH^T,) is ^(Ei<i^ , < m ' l < ri / )• 

Theorem 10. Fix some odd d > 3. There exists a set S of spheres in ~E d , consisting of rii spheres 
of radius pi, 1 < i < m, with p\ < P2 < ■ ■ ■ < p m an d m > 3 fixed, such that the complexity of the 

L-J 

convex hull CH d (E) is ^(J2i<i^j< m n i n j 2 )• 

Consider again the injection cp : CHdiJ^) —> V . We have shown above that the worst-case 
complexity of C^(S) is ^(X^i<i^<m n i n j 2 )) when d > 3 is odd and m > 2 is fixed. Since ip is 
injective, this lower bound also applies to the complexity of V . This establishes our lower bound 
claim in Theorem [6j 

Corollary 11. Let & = \V\ , V2, ■ ■ ■ , V m } be a set ofm d-polytopes, lying on m parallel hyperplanes 
of~E d+l , with d > 3, d odd, and both d and m are fixed. The worst-case complexity of CHd+i(&) 

ts tt{J2i<i^j< m n i n \ 2i ), where rii = fo{Vi), 1 < i < m. 

2 Combinatorial equivalence here means that each face of circularity £ of A r corresponds to a unique (d — £~ l)-face 
of A . 
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5 Computing convex hulls of spheres 



In this section we focus our attention to the computation of the convex hull CH^E) of E. We use 
the same notation as in Section [JJ 

Polarity with respect to the point O is a one-to-one transformation which maps points of E d+1 
distinct from O' to hyperplanes of E rf+1 that do not contain O' (e.g., see [8j Section 7.1.3]). For a 
point q in E rf+1 distinct from O', we define its polar hyperplane q* by 

?* = {!£ E d+1 \q-X = 1}, 

while for a hyperplane H of K d+1 not containing O', the polar point or po/e H* of is defined by 

H* ■ X = 1, \/X £ H. 

The po/ar sei of a set of hyperplanes is the set of the poles of these hyperplanes. Finally, we denote 
by H~ the closed halfspace bounded by H and containing O'. 

Consider again the point set P in E rf+1 that we get by lifting the spheres of S. The polytope 
V* in K d+1 defined as the intersection of the halfspaces Pk*~ , 1 < A; < ra, is the polar polytope or 
dual polytope of V . V and V* are dual in the sense that there is a bijection between the I- faces of P 
and the (d — £)-faces of V*, which is inclusion-reversing: a hyperplane supporting V along an £-face 
F, has its pole on the corresponding (d — ^)-face F* of V*. In a similar manner, we can define the 
polytope P*, which is the polar polytope of V . 

Let Ho' be the hyperplane parallel to Hq passing through O' and let Hq, denote the open 
halfspace of K d+1 above Hq 1 - Boissonnat et al. O Proposition] have shown that: 

1. The polar set of the hyperplanes, which are translated copies of the hyperplanes tangent to 
Ao, is the hypercone K with apex at O 1 , vertical axis, and angle at the apex equal to \. 

2. The polar set of the hyperplanes above O' is the halfspace Hq,. 

An immediate consequence of these properties is that the polar set of the hyperplanes supporting 
V, tangent to at least one hypercone of A along a generatrix, and above O', is V* fl^Tl Hq,. 

Boissonnat et al. [6] have used this property in order to propose an algorithm for computing 
C-ffrf(S) in 0(n^2^ + nlogn) time for any d > 2. Below, we describe a slightly modified algorithm 
that takes into account the fact that the radii of the spheres in £ can take a fixed number of m > 2 
distinct values. Our algorithm consists of the following five steps: 

1. For all i with 1 < i < m: determine the set Pi = P Pi Ilj, construct the convex hull V% = 
CHd(Pi), and compute the set Pi. 

2. Compute the polytope V = CHd+i(P), and choose a point O' inside P. 

3. Compute the polar polytope V* of V with respect to O' . 

4. Compute the intersection between P*, the hypercone K, and the halfspace Hq,. 

5. Compute the incidence graph of the facets in CH^iT,) from the incidence graph of the faces 
of P* intersecting K and Hq, . 

L-J 

Determining all the sets Pi takes G(n) time, whereas constructing the polytope V% takes 0(n i 2 + 
rii log raj) time. Determining Pi from from the representation of Vi can easily be done in 0{rii log raj) 
time. We thus conclude that step 1 of the algorithm takes 0(n^i + nlogn) time. Computing V* 
from P takes linear time in the size of P. Moreover, the complexity of computing P* flKfl Hq, 
is linear in the complexity of P*, and thus P, since K and Hq, have constant complexity. The 
complexity for computing the incidence graph of CHd(T,) from the incidence graph of P* n K C\ Hq, 
is again linear in V* and V . Therefore, the total time needed for steps 2-5 of the algorithm is a 
linear function of the time needed to compute V . Summarizing: 
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Theorem 12. Let J] be a set of n spheres in E , having m distinct radii pi, p%, . . . , p m , with 
d > 3, d odd, and both d and m are fixed. Let m be the number of spheres in S with radius pi, 
1 < i < m. We can compute the convex hull C-ffrf(S) in 0(nW + nlogra + Td+i(ni,n2 ■ ■ ■ i n m)) 
time, where Irf+i (m , ri2 ■ ■ ■ , n m ) stands for the time to compute the convex hull of m d-polytopes 
{V\,V2, ■ ■ ■ ,V m } lying on m parallel hyperplanes in K d+1 , with m = fo(Vi), 1 < i < m. 

As described in Section [3l 

L~J 

T d+l {m, . . . ,n m ) = 0{(J2i<i^j<m n i n j 2 
for any odd d > 5, whereas for d = 3 we have 

r 4 (m, ...,n m ) = 0(min{n 4 / 3+e + (Ei^< m ^nf ) log n, (£i<.^< m ^ J ) ^g 2 n}). 

I -J 

Hence, we can compute the convex hull Ci?d(S) in 0((Y2i<ijLj< m n i n j 2 )l°g n ) time for any odd 
d > 5, and in 0(min{n 4/3+e + (X)i<#j< m «i^- 2 ) log n, (Z)i<i^j<m n * n j 2 ) lo g 2 n i) time for ^ = 3 - 

6 Summary and open problems 

In this paper we have considered the problem of computing the worst-case complexity of the convex 
hull CH c [ + i(^ > ) of a set & = {V\, V2, ■ ■ ■ ,V m } of m convex d-polytopes lying on m parallel hyper- 
planes of E d+1 , for any odd d > 3. Denoting by rii the number of vertices of Vi, we have shown 

I -J 

that the worst-case complexity of CHd+i(3 g ) is 0{^2 ll< -_ L - <m nin- 2 ). This result suggests that, in 
order to compute CH ( i+i{3 ) 1 it pays off to apply an output-sensitive convex hull algorithm to the 
set of vertices in 2? . Indeed, for any odd d > 5 Seidel's shelling algorithm |27| . or its modification 

by Matousek and Schwarzkopf |21] . results in a 0({^2 1< j_ L j <m nin- 2 )logn) time algorithm, where 
n = Yli=i n i- F° r d = 3. the divide-and-conquer algorithm by Chan, Snoeyink and Yap |11| can 
be competitive against Matousek and Schwarzkopf's algorithm; hence, we may compute CH^{&) 
in 0(min{n 4/3+e + (Ei<^j< m ?VV) logn, (Ei<^< m njnj) log 2 n}) time, for any fixed e > 0. The 
above algorithms are nearly optimal for any odd d > 3; it remains an open problem to compute 

CHd-\-i{&) in worst-case optimal 0(Y^i<j^j<m n i n j 2 + ?^logn) time. 

A direct consequence of our bound on the complexity of CHfi + i{^) is a tight asymptotic 
bound on the worst-case complexity of the (weighted) Minkowski sum of two polytopes in any 
odd dimension d > 3. More precisely, consider a n-vertex (f-polytope V and a m-vertex d-polytope 
Q, and embed them on the hyperplanes {x^+i = 0} and {x^+i = 1} of E d+1 , respectively. The 
weighted Minkowski sum (1 — X)V © XQ, A G (0, 1), is combinatorially equivalent to the intersection 
of CHd+i^V, Q}) with the hyperplane {xj+i = A}, whereas the Minkowski sum V © Q is nothing 
but © |Q, scaled by a factor of 2. Applying our results, we deduce that the complexity of 

(1 - X)V © AQ (resp., V © Q) is e(nm^ + 

mnlaJ) for any odd d > 3. We would like to extend 
this tight bound to (weighted) Minkowski sums where the number of summands is greater than 2. 

Capitalizing on our result on the complexity of convex hulls of convex polytopes lying on parallel 
hyperplanes, we have shown that the worst-case complexity of the convex hull CH^CE) of a set E 

of n spheres in K d with a fixed number of m distinct radii pi,p2, ■ ■ ■ ,p m I s ^(X^Kj^j<m n i n j 2 )> 
for any odd d > 3, where rtj is the number of spheres with radius pi. By means of an appropriate 
construction, described in Subsections 14.11 and 14. 2\ we have shown that the upper bound above 
is asymptotically tight, implying that our upper bound for CH ( i + i(3 g ) is also tight. By slightly, 
but crucially, modifying the algorithm of Boissonnat et al. [6], CHd(T,) may be computed in 
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0(n*-2i + nlogn + Td+i(ni, . . . ,n m )) time, where 1^+1(711, . . . ,n m ) stands for the time needed to 
compute the convex hull of m d-polytopes lying on m parallel hyperplanes in E rf+1 , where the i-th 
polytope has rtj vertices (cf. Section [5]). 

Finally, Boissonnat and Karavelas [7] have shown that convex hulls of spheres in K d and addi- 
tively weighted Voronoi cells in K d are combinatorially equivalent. This equivalence suggests that 
we should be able to refine the worst-case complexity of an additively weighted Voronoi cell in any 
odd dimension, when the number of distinct radii of the spheres involved is considered fixed. 
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